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Abstract. A Monte Carlo simulation of mergers in clus- 
ters of galaxies is carried out. An "explosive" character of 
the merging process (an analog of pha se transition), sug- 
gested earlier by Cavaliere et al. ( 1991 ), Kontorovich et al. 
(1992), is confirmed. In particular, a giant object similar 
to cD-galaxy is formed in a comparatively short time as a 
result of mergers. Mass and angular momentum distribu- 
tion function for galaxies is calculated. An intermediate 
asymptotics of the mass function is close to a power law 
with the exponent a « 2. It may correspond to recent 
observational data for steep faint end of the luminosity 
function. The angular momentum distribution formed by 
mergers is close to Gaussian, the rms dimensionless angu- 
lar momentum S/(GM 3 R) 1 ^ 2 being approximately inde- 
pendent of mass, which is in accordance with observational 
data. 



Key words: galaxies: clusters: general - galaxies: cD - 
galaxies: interactions - galaxies: mass function - galaxies: 
statistics 



1. Introduction 

In recent years, evolution of galaxy clusters has attracted 
ever more attention. On the one hand, it has become acces- 
sible for observations, both with ground-based and space 
instruments. On the other hand, discovery of Butcher- 
Oemler effect testifies to an epoch of fast evolution in 
clusters.- associated with galaxy interaction. Mergers of 
galaxieso are considered to be one of the most probable 
explanations for the change of colour which accompanies 
appearance of ellipticals and lcnticulars instead of e arly 
type spirals at z ~ 0.2-0.4 (see, e.g., Dressier et al. 1994 



and cited there) . Comparatively fast evolution of clusters 
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1 Another possible ex plana tion is "galaxy harassment" dis- 



and groups, caused by mergers, is also confirmed by the 
results of direct numerical simulation (see, e.g., Barnes 



1996|) 



One of the effects associated with mergers is rapid evo- 
lution of the galaxy mass and angular momentum distri- 
bution function f(M, S, t), related to appearance of mas- 
sive galaxies. Calculation of f(M, S, t) is of great interest 
both by itself and in the context of a merger model of 
activity suggested by Kats & Kontorovich (1990, 1991). 
Given f(M, S,t), one can predict the luminosity function 
of active galactic nuclei in this model. 

The mass function of merging galaxies $(M, t) can be 
described in terms of the Smoluchowski kinetic equation. 
As Cavaliere et al. (1991 1992) and, independently, Kon- 



torovich et al. ( 1992 ), Kats et al. ( 1992 ) have shown, merg- 
ers result in an analog of phase transition (or "explosive 
evolution"), which implies fast formation of a distribu- 
tion tail (corresponding to massive galaxies) and a "new 
phase": cD-galaxiesc. 

The joint distribution function which takes into con- 
sideration both mass and angular momentum can be de- 
scribed by a generalized Smoluchowski equation (see Kats 
& Kontorovich 199C , 1992 , where this equation was solved 
in the simplest case of a constant merger probability - 
an analog of phase transition does not take place in this 
variant - and without allowance for the orbital angular 
momentum). In this paper we present the results of sim- 
ulation of galaxy mass and momentum evolution in clus- 
ters, with the orbital momentum and more realistic mass 
dependence of the merger probability taken into account. 

In Sect. Hwe discuss Monte Carlo simulation of merg- 
ers. In SectHs we compare the results for $(M, t) with a 
direct numerical solution of the Smoluchowski equation. 
Section H contains discussion of the results. 



This transition has been known for a long time in oth er ap- 



cussed by Moore et al. (1996) 



plications of the Smoluchows ki equ ation ( see St ockmayer 
Trubnikov |l97l| ; Voloshchuk [L984| ; Ernst |l986| ). 
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2. Joint mass and angular momentum distribution: 
Monte Carlo simulation 

2.1. Generalized Smoluchowski equation and merger prob- 
ability 

The kinetic equation which describes f(M, S, t) (the gen- 
eralized Smoluchowski equation) is 

§£ = J W MSlMlSlM2S J 1 f 2 dM 1 dM 2 d 3 S 1 d 3 S 2 

-2 j WM 1 s 1 \MSM 2 sJf2dM 1 dM 2 d 3 S 1 d 3 S 2 , (1) 

where f x = f{Mx,S x ,t) etc. The function W M s\m 1 S 1 ,m 2 S 2 
(the kernel of the equation) is a characteristic of the prob- 
ability for the merger M 1 Si,M 2 S 2 -> MS. Taking into 
account mass and momentum conservation laws, W can 
be rewritten as 



W, 



MS\M 1 S 1 ,M 2 S 2 



= u 



MS\M 1 S 1 ,M 2 S 2 



x8(M - Mi - M 2 )S(S -Si- S 2 ) 



(2) 



(without orbital momentum) . If we take into consideration 
the orbital momentum J, the second ^-function should be 
replaced by a function of a finite width ~ J. The function 
U can be calculated as U = av, where a is the merger 
cross-section, v is the relative velocity at infinity, the bar 
means an average over velocities. 

An exact computation of collision dynamics and de- 
termining the merger cross-section is a very complicated 
problem. Nevertheless, main features of this process are 
known, both from analytical consideration and numeri- 
cal experiments (Roos & Norman 197S| ; Aarseth & Fall 



1980; Farouki & Shapiro 1982; Farouki et al. 1983; Chat 



terjee 1992) and enables one to formulate conditions nec- 
essary for a merger: the galaxies must pass at a small 
distance (interaction is especially intense if the outer 
parts overlap) and the relative velocity must be small 
enough. Namely, we shall assume below that a merger oc- 
curs provided that (i) the minimal distance between the 
two galaxies is less than the sum of their radii Ri + R 2 
and (ii) the relative velocity at infinity is less than some 
limit value which is of the order of the escape velocity 
v e = y /2G(M 1 +M 2 )/(Ri+R2), i-e., v < C«g, C ~ !• 
Taking into account gravitational focusing, we may de- 
rive from the former condition that the impact parame- 
ter poo < (Ri + R 2 )^J\ + v 2 /v 2 . Thus, the merger cross- 
section is 



^(i?!+i? 2 ) 2 (l+^A' 2 ), V<(v g 



0, 



v > (v g . 



(3) 



We shall assume that the galaxy peculiar velocity distri- 
bution is Gaussiand with the root mean square velocity 



3 It is essential for the further consideration that t> rms does 
not depend on mass. Such a behaviour is typical for gravitating 
systems which have passed through the violent relaxation stage 
(Saslaw [L987|). 



u rms (obviously, the relative velocity distribution (at in- 
finity) in this case is also Gaussian, but the root mean 
square velocity is y/2v ims ). So, 



fC v e 

U = av= n(Ri + R 2 ) 2 {1 + v 2 J v 2 )v 
Jo 

x (47T/3)- 3 / V m 3 s exp(-3 V 74 Ur 2 ms ) W dv. 
After integration we obtain: 
(7 = 4( 7 r/3) 1 / 2 z; rms (i? 1 +i? 2 ) 2 



(4) 



A + l 



-C 2 a 



(( 2 A + A+1) 



A 



(5) 



If the density of a galaxy p ga \ does not depend on 
its mass, then the radius R and mass M are related as 
R oc M 1 / 3 . For Faber-Jackson and Tully-Fisher laws 
(L oc V 1 ), using the virial theorem (MV 2 cx GM 2 /R) 
and the mass-luminosity relation of the form L oc M , we 
obtain R oc M 1 / 2 . Below we shall take 



R = CM?, where (3 = 1 - 1 



3 2' 



(6) 



U 



Asymptotical behaviour of Eq. (5) is 

c 2 (Mi+M 2 ) 2 , M<M b 
ci +/9 (Mi + M 2 )(Mf + Mi), M > M b , 



(7) 



the coefficients c 2 = 3(3tt) 1 /2((2 + ^/2)G 2 /v^ ms , c 1+/3 
2(37r) 1 / 2 CG/w rms . Here mass Mb corresponds to v rmB 



M b - (Cw r 2 ms /G) 1/(1 -' 3) - 10 9 -10 



10 2/(1-/3) 



M Q , 



(8) 



where v 7 = w rms /10 7 cm s _1 . We shall assume C = 1 in 
the further consideration. 

Note that the difference of U in the small mass re- 



gion from that in Cavaliere et al. 1991 is caused by the 
velocity restriction for mergers in Eq. (g). For M < M b 
collisions without mergers are more probable than mergers 
(a = tt(Ri + R 2 ) 2 )- The latter, however, gives the "explo- 
sive" evolution in the small mass region too. Effects re- 
lated to collisions without mergers may be important (see 
footnote |l|), however, this question is beyond the scope of 
this paper. 

As v Tms is, in general, a function of time, the coef- 
ficients c 2 and C1-1-/3 and m ass M b depends on time too 
(see, e.g., Kontorovich et al. 1992 ). Below we shall assume 
them to be constant, neglecting changes of velocities and 
masses both due to capture of new members or evapo- 
ration of galaxies from clusters and groups, and mergers 
itself, or other reasons. 

In the region M ^> Mb galaxy peculiar velocities are 
much less than stellar velocities in the galaxy; the rela- 
tionship for M -C Mb is inverse. In this paper we consider 
both the asymptotical regions M 3> M b and M <C Md 
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and the intermediate case where Mb is within the range 
of simulated masses. 

Numerical experiments of galaxy merging (Farouki & 
Shapiro 1982| ) show that the merger probability U de- 
pends both on masses and momenta, reaching a maximum 
when Si, S2 and J have the same direction. Nevertheless, 
this dependence is less essential than the dependence on 
masses (oc M", u = 1 + /3 for M » M b , which leads 
to the "explosive" evolution). In this paper we shall use 
the simple model (||) which does not take into account 
the dependence of the merger probability on the mutual 
orientation of angular momenta. 

It is known from the Smoluchowski equation theory 
(Ernst |1986| ; Voloshchuk |l984| ) that behaviour of the solu- 
tion essentially depends both on the homogeneity power 
u in the mass dependence of U and on the asymptotics of 
U for very different masses which can be characterized by 
exponents U\ 2- 

U{M ll M 2 ) « c u M^M^ , Mr < M 2 , u x + u 2 = u.(9) 

For the case of galaxy mergers (Eq. (^)), obviously, u\ = 
and u 2 = u = 2 (Af < M h ) or 1 + (3 (M > M h ). If u > 1 
then an analog of phase transition takes place in the sys- 
tem. For an initial distribution localized in the small mass 
region, a slowly decreasing distribution tail is formed; in a 
finite time the tail reaches infinite masses (in an idealized 
case). The second moment of the distribution becomes 
infinite at t — t CI . This phenomenon may, in principle, 
lead to fast formation of m assive galaxies and quasars by 
mergers (Kontorovich et al. 1992 ), and also to forma tion o f 
cD-galaxies in groups and clusters (Cavaliere et al. 1991 ). 



2.2. Numerical simulation of mergers: methods 

Direct analytical or numerical solving the generalized 
Smoluchowski equation with orbital momentum is a very 
complicated problem due to the great number of variables. 
This difficulty can be avoided by using numerical Monte 
Carlo simulation of merging process. In this section we 
present such simulation. A finite system consisting of A 
galaxies (referred to as "cluster" below, though it may also 
be a group) was considered. Pairs of these galaxies merged 
(with probability, proportional to U(M 1 ,M 2 , Si,S 2 )) un- 
til the number of the galaxies reduced to some Af . For 
each merger mass and angular momentum conservation 
laws M = Mi + M 2 , S = Si + S 2 + J were fulfilled. 
Distribution of the angular momenta (intrinsic Si i2 and 
orbital J) over directions was taken isotropic. It was as- 
sumed that the merger probability U depends only on 
masses according to Eqs. (5) and (|^) and does not de- 
pend on momenta. The absolute value of the orbital 
momentum was computed in accordance with Eq. (|J), 
namely, J = m^m 2 v Poa , v being a random number 
distributed in the range to v g with probability den- 
sity f(v) oc v 2 exp(— 3v 2 / 4v^ ms ) , and the impact parame- 
ter poo being a random number distributed in the range 



to (i?i + R 2 )(l + vl/v 2 ) 1/2 with probability density 
fipoo) oc poo- The initial galaxy mass distribution was 
chosen exponential $o(M) oc e~ M ' M ° (the results are in- 
dependent of the exact expression for initial distribution, 
only the fact that it is localized in the region of small 
masses ~ Mo and decreases rapidly for large ones is es- 
sential). Instead of the momentum S, it is convenient to 
introduce a dimensionless momentum A = MR ( 2 G^i/Ry/ 2 ' 
similar to Peebles' parameter (SE 1 / 2 )/(GM 5 / 2 ) ( see also 



Doroshkevich 1967). Dimensionless momenta A of the ini- 
tial galaxies were distributed uniformly in the range to 1 
in our simulation (as for <&o(M), the exact form of the dis- 
tribution is unessential). Computations were carried out 
for the following parameters: (3 = 3—5, A = 10 3 -10 5 , 
A f = (l-Tl-lO^A. 

In conclusion we describe the procedure of simulation. 

1. A initial galaxies are simulated, with mass and angular 
momentum distributed according to fo(M, S). 

2. Two random integer numbers, i and j, distributed uni- 
formly in the range [L,n] (n is the current number of 
galaxies) and satisfying the condition i =/= j, arc simu- 
lated. 

3. Galaxies number i and j merge with probability p = 
U(Mi, Mj)/U max . The mass and momentum of the 
new galaxy are calculated as described above. With 
probability (1— p) the galaxies do not merge, and jump 
to item 2 is executed. Here U max = max U(Mi, Mj) 

l<ij<n 

depends on time. 

4. Items 2-3 are executed until the number of galaxies n 
becomes equal Af. 

The algorithm used in our simulation was somewhat 
different from the simplified scheme given above. The rea- 
son was that the merger probability p = U (Mi, Mj)/U max 
is very small for the majority of galaxies and so the simula- 
tion time is very large. In actual simulation the algorithm 
was modified as follows: 

— the probability to choose the z-th galaxy in item 2 was 
ki/N instead of 1/A, where fcj is the number of initial 
galaxies which have subsequently merged into the i-th 
galaxy; 

— to compensate this change, the function U'(Mj,, Mj) = 
U(Mi,Mj)/(kikj) was used instead of U(Mi,Mj) in 
item 3. 

Obviously, these modifications do not influence the result 
of the simulation. In the same time, the number of cycles 
reduces because the average value of U' /t/j' nax is closer to 1 
than the average value of U/U max . 

2.3. Numerical simulation of mergers: results 

After some time, an analog of phase transition takes place 
in the system of merging galaxies, similarly to what oc- 
curred in the work by Cavaliere et al. (1991). The system 
divides into two phases: a giant galaxy which contains a 
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~~ i 1 r~ 

U oc (Mi + M 2 ) 2 



logio A 




login M 

Fig. 2. Masses M and dimensionless angular momenta A 
of simulated galaxies for the case M <C Mb. N = 10 4 , 
N{ = 4-10 _1 A^, 100 clusters (4-10 5 galaxies). 



major part of the total mass and many small galaxies (in 
the cases u = 3/2 and u — 2 this transition is more evident 
than for u = 4/3). These giant galaxies formed by mergers 
can probably be identified as cD-galaxics in the centres of 
groups or clusters (cf. Hausman & Ostriker 1978; Cava- 



liere et al. 1991). The majority of low-mass galaxies are 
those which have never merged (Fig. [j]). Such behaviour 
is related to strong dependence of the merger probability 
on masses (u — 1 + (3 or u — 2) which, in turn, results in a 
steep mass function (a ~ 2 for u — 1 + /3, a ^ 2 for u = 2, 
see below). 



H 1 1 1 r~ 

U oc (Mi + M 2 )(M 1 1/2 + Ml' 2 ) 




logic M 

Fig. 1. Masses M and dimensionless angular momenta A of 
simulated galaxies (M ^> Mb): a initial, b and c formed by 
mergers (each dot represents one galaxy). A distribution tail, 
independent of the initial conditions, is formed due to mergers. 
The right-hand part of the figure (M ~ 10 4 ) corresponds to 
cD-galaxies. For u = 3/2 (Fig. c) separation of galaxies into 
the two phases, normal and cD, can be seen better than for 
u = 4/3 (Fig. b). Each diagram shows 10 5 galaxies (b and 
c - 100 clusters of N { = 1000 galaxies each); N f = 10 _1 7V, 
TV = 10 4 . 



login* 




log xo M 

Fig. 3. Mass function obtained by Monte Carlo simulation 
(M 3> Mb). The values of parameters are the same as in Fig. |lj 
Mass is given in units Mo, $(M) is normalized to unity. The 
part of the plot near M = 10 4 corresponds to cD-galaxies. 



The mass function formed by mergers in the case U oc 
(Mi + M 2 )(Mf + Mjf ) (i.e., M > M b ) is shown in Fig. § 
In the region M <C M <C M max it is close to a power law 
M~ a , a s» 2. The rise near the right-hand end corresponds 
to cD-galaxies the masses of which are comparable to the 
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distribution 
-obtained 

distribution 



logic Af 



Fig. 4. Mass function obtained by Monte Carlo simulation for 
different M b (N{ = lO^Af). The average slope changes with 
M b . 



total masses of their clusters. In the case U cx (Mi + M2) 2 
(i.e., M <C Mb), and in the intermediate case of a finite 
Mb the mass function is steeper. 

The obtained momentum distribution at fixed mass 
in the asymptotical region of large masses is close to the 
normal distribution (Fig. |^a). Thus, the distribution tail 
may be represented as 



0.5 I 1.5 2 2.5 

^2 S 



3 Sri 



N(< S) 

iVtot 



/(M, S) » 3(M) 



27 V 

T 



-3/2 



6X P _ o^T 



3S 2 
2^2 



(10) 




where <£>(M) is the mass function, S 2 (M) is the average 
square value of the momentum 5* for given mass M. cD- 
galaxies in the cases u = 3/2 and u — 2 make an exception 
(Fig. ^b). Kolmogorov-Smirnov test shows evidently that 
their momentum distribution differs from the normal one. 



2 S 



3 St: 



Fig. 5. Comparison of the angular momentum cumulative 
distribution at a fixed mass with the normal distribution 
(Eq. ([[(})) for M > Mb- a The region of the distribution tail 



For U cx (Mi + M 2 ) (Mf + M|) simulation shows that, ^10 M = 1.8) for u = 3/2. The distribution is close to the 



irrespectively of the initial momentum, the root mean 
square value of the dimensionless momentum A rms be- 
comes constant at large masses: 



A,. 



const ~ 0.1, 



that is 



S* 2 (M) ~ 0.01 



MR 



2GM 
R 



(11) 



(12) 



(Figs, ^a and |T]). The distribution at small masses depends 
on/ (M,S). 

For U oc (Mi + M2) 2 the dimensionless momentum 
decreases with mass (Figs. ||b and ^|). 

The fact that A rms « const for U cx (Mi + M 2 )(Mf + 
Mf) at large masses has a simple qualitative explanation. 
Consider the change of the mass and momentum due to 
mergers. As the mass function decreases rapidly (i.e., the 
number of small galaxies is very large) and u\ = 0, it 
is natural to suppose that the main contribution to the 



normal one. b cD-galaxies for u = 3/2. The distribution dif- 
fers from the normal one (the significance level in the Kol- 
mogorov-Smirnov test is ~ 10 -9 ). 



change of the mass M and momentum S of a given massive 
galaxy is associated with accretion of small galaxies ~ 
Mq <C M (as, e.g., in Kontorovich et al. 1992 ). However, 
in this case the situation is different. The rate of changing 
M and S due to mergers with low mass galaxies (< M) 
can be expressed as 



M 



M 



U(M, Mi)Mi$(Mi) dMi 



cx M 



M 

"- I M 111+1 




$(Mi)dMi, 



(13) 



M 



S 2 = / U(M, Mi)(S' 2 (M 1 ) + J 2 (M, Mi))$(Mi) dM x 
Jo 

r u i+ 2 tf>fi\/T^\AH/r, (14) 

'0 



xM u 2 +i+0 1 M 1 tl+:i $(Mi)dMi, 
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Fig. 6. Dependence of the root mean square dimensionless 
momentum A rms on mass, a M 3> Mb. In the left-hand part 
of the plot the main contribution is given by small galaxies 
which have never merged, A rmB being determined by the initial 
distribution for them. For large galaxies momenta are deter- 
mined by mergers and do not depend on the initial conditions, 
Arms ~ const, b The same for different M b (JV f = 10 _1 iV). 
A rms decreases with M. 



since S 2 {M 1 ) < J 2 (M,Mi), J 2 (M,Mi) ~ M 2 GMR for 
Mi <C M. In our case u\ — 0, the slope of the mass 
function a ~ 2. Therefore, the main contribution to inte- 
gral (14) is given by large galaxies (~ M); contribution of 
large and small galaxies in Eq. (13) are of the same order. 

Next we consider successive mergers of galaxies with 
equal large masses M and momenta S ims (M). n The new 
galaxy has a mass M' = 2M and momentumQ (S 1 ) 2 = 
2S 2 + J 2 = 2S 2 + MGMR. Thus, the new value of the 

2 -2-/3 A 2 + 2-4-/3. Afj 

ls tends to the equilib- 



dimensionless momentum (A') 2 = 
a result of successive mergers, A r 



rium value A r 



i!" 2 4 - 2 % w - 23 ~ - 25 (/? = H)- 



So, the result A 2 ss const is a consequence of the fact 
that the main contribution to the change of mass and mo- 



4 The orbital momentum J 2 = (MM/(M + M)) 2 v 2 pL, the 
impact parameter pJL = \{R+ R) 2 v1/v 2 , so, J 2 = MGMR. 



mcntum is given by mergers between comparable mass 
galaxies. The obtained value for A rms is somewhat differ- 
ent from Eq. (|ll]) due to simplifying assumptions made 
in the derivation. Note, that for the isotropic momentum 
distribution without allowance for the orbital momentum 
S 2 oc M, A 2 oc M~ 2 ~^; for the anisotropic distribution 
(Kats & Kontorovich |l990i [L99lD S oc M, A oc M-( 1+/3 )/ 2 . 

In the case J7 oc (Mi + M2) 2 the main contribu- 
tion to Eqs. (13), (14) is given by small masses ~ Mo 
(due to large a) and we can expect decreasing of A rms 
MS M" 1 / 2 . This can be demonstrated as follows. If inte- 
gral J Q M" 1+2 $(M)dM converge at infinity then it is 
possible to replace the upper limit of integrals (13), (14) 
by infinity. Then 



A/ocM" 2 , S 2 oc M U2+1+f> . 
Therefore, 



(15) 



dS 2 
dM 



M 1+ P S 2 oc M 2+f \ A 2 



,S* 2 



M 3 +/3 



oc M -1 . (16) 



An analysis of observational data (Kontorovich & 
Khodjachikh 1993 Kontorovich et al. [1995a ) confirms 
that SVms oc M fe , the coefficient k being rather close to 
(3 + /3)/2, which is in accordance with Eq. Jl^). 

Note that allowance for dependence of the merger 
probability on momenta may give an increase of A 2 : nu- 
merical experiments show that a merger is more probable 
if all momenta have the same direction (Farouki & Shapiro 



1982 ). On the other hand, A 2 is sensitive to the exact form 
of the merger cross-section (in particular, to the value of 
C in Eq. (|). 

3. Comparison of simulation results with solution 
of the Smoluchowski equation 

Integrating the generalized Smoluchowski equation (1) 
over momenta, one can obtain the ordinary Smoluchowski 
equation which describe the evolution of the galaxy mass 
function: 

d$(M, i) 

Ft 

M 

C/(Mi, M - Mi, t)$(Mi, t)$(M - Mi, t) dMi 

r OC 

-2$(M,f) / [/(Mi,M,i)$(Mi,f)dMi. (17) 
J 

Solving this equation is another independent way to 
find <&(M, t). In this section we compare the results ob- 
tained by Monte Carlo si mulation with the obtained ear- 
lier (see Kontorovich et al. 1995b and Krivitsky 1995 ) re- 
sults of direct numerical solving the Smoluchowski equa- 
tion with kernels U oc (Mi + M 2 )(Mf + M|) and U oc 
(Mi + M 2 ) 2 . 

For numerical solution of Eq. (17) and analysis of the 
obtained results we used methods described in Krivitsky 
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galaxies, obtained from the solution of the Smoluchowski 
equation, is shown in Fig. |^. The moment displayed in 
Fig. |l|b (JVf = O.lJVo) approximately corresponds to t « 
0.3 (whereas t CT w 0.26). 

In the case U = c 2 {M 1 + M 2 ) 2 (Fig. 0a) the interme- 
diate asymptotics is not as close to a power law as for 
U = c u (Mi + M 2 )(Mf + Mjf ). An effective slope a in this 
case is 2-3. The value of t CT is ~ 0.02, that is the phase 
transition is very fast. However, in this case the time de- 
pendence of the distribution moments is non-power and 
t cr cannot be determined accurately (Kontorovich et al. 



19951:; Krivitsky 1995) 



N(t) 
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0.2 I- {Ml 1 3 + M 1/3 



U oc (Mi + M 2 )x 
2 1 



J_ 
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tcr 0.3 



logi M 



Fig. 7. a Numerical solution of the Smoluchowski equation 
in the small mass region (M max = 10 4 , M is measured in 
units M , time in units l/{c u N M%), $ in units N /M ). 
b Solution in the large mas s regi on, for comparison (the same 
case as in Cavaliere et al. (1992), but in a wider mass range 
and somewhat more precise). 



Fig. 8. Time dependence of the number of galaxies, obtained 
from the solution of the Smoluchowski equation (upper time 
scale for it = 2, lower time scale for u = 4/3). Though the total 
number of galaxies decreases, at any t the number of low-mass 
galaxies which have never merged is much less than the number 
of massive tail galaxies. 



1995; Kontorovich et al. 1995b. Figure ^b shows a plot 
of the obtained mass function for U = c u {Mi + M2) x 
(Mf +M% ). One can see that an intermediate asymptotics, 
close to a power law M~ a , is formed in the region Mo 
M <C M max . The exponent a w 1.9 for u = 4/3 and 
a w 2.1 for u = 3/2 (in the latter case a is defined worse, 
see below). The value 1.9 differs from a ~ 2.15 given by 
Cavaliere et al. Jl992 ) but, in our opinion, agrees rather 
well with the .plot in their work. The value of t cr for the 
two cases areQ, respectively, 0.26 and 0.1 (in the same 
units as in Fig. 0) for the initial distribution $>q{M) = 
{No/Mo)e~ M / M ° . The time dependence of the number of 



5 Since a characteristic time associated with mergers is of the 
order of {cuNoMq) -1 , such a value of t CT means that the phase 
transition is comparatively fast. Due to considerable contri- 
bution of a comparatively small number of appearing massive 
galaxies, the time t cr corresponding to the "phase transition" 
is much less than a characteristic time (trim) -1 , where a is the 
merger cross-section for typical galaxies, n is the concentration 
of such galaxies, v is the average velocity. 



In numerical solving Eq. (17), a finite limit mass M max 
was introduced: the integral from to 00 in the right-hand 
part was replaced by the integral from to M max . Phys- 
ically, such a substitution corresponds to a sink at large 
mass es. As i t was shown by Krivitsky 1995 , Kontorovich 
et al. |l995b| , consequences of this replacement are differ- 
ent for kernels with 112 < 1 and u-i > 1. In the case of 
kernel ((7J) which belongs to the class Ui > 1, existence 
of M max and its value essentially influence the solution, 
in particular, the number of galaxies N — J Q $ dM and 
distribution moments M. ^ = J °° $M P dM as functions 
of time, the value of t CI etc. Moreover, van Dongen (1987) 
showed that the limit M max — > 00 does not exist at all 
for 112 > Ui + 1 (this is the case for Eq. ([?])). The influ- 
ence of M max increases as ui becomes farther from 1 : for 
U tx (Mi + M 2 )(Mf + Mf), especially if (3 = 1/3, this 
influence is moderate. The farther ui moves from 1, the 
greater the difference between the intermediate asymp- 
totics and the power law becomes and the worse defined 
t CI and a are; this is the case for U cx (Mi + M2) 2 . 
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As known (see van Dongen & Ernst 1988; Voloshchuk 



1984 ; Krivitsky |1995[ Kontorovich et al. |l995b| ), in many 
cases the solution for U2 < 1 is self-similar: $(Af) w 
p~ T (t)ip(M / p(t)) for M ^> Mo, where <p is a time- 
independent function, describes a "front" moving to 
greater masses. The numerical solution shows that $(M, t) 
for U oc (Mi + M 2 )(A/f + Aff ) is closer to the self-similar 
form for lower f3. The mass function for U oc (Mi + M2) 2 
is not self-similar (nonlocal case, see discussion in Kon- 
torovich et al. 1995b| ). However, for t > t CI the shape of the 
curve becomes nearly constant (Fig. 0a): a different self- 
similarity appears, because mergers with the cD-galaxy 
dominate and the dependence of U on the smaller mass 
vanishes. 

Both in direct solution of the Smoluchowski equation 
and in numerical simulation of mergers, there exists a fi- 
nite limit mass: the mass of the sink in the former case, the 
total mass of the system in the latter one. However, the 
problems which are solved in this section and in Sect. [2] 
are not equivalent. Nevertheless, simulation shows that, in 
spite of the essential influence of the limit mass (112 > 1), 
the mass function obtained by simulation of mergers in 
Sect. U and the one obtained as a direct solution of the 
Smoluchowski equation have good agreement in the region 
M <C M max . So we can make a conjecture that <!>(M) for 
M <C M max depends only on the value of the limit mass 
and does not depend on its nature. 

4. Discussion 

Simulation confirms the possibility that massive galaxies 
may form by mergers, moreover, this process has an "ex- 
plosive" character and is an analog of phase transition, 
cD-galaxies (with mass comparable to the total galaxy 
mass of the cluster) being formed as a new phase. What 
are the conditions which make this process possible? The 
expression for t cr may be written as 

t C r = — S . ,,. (18) 
c u N M^ 



galaxies to the average density of the Universe; the co- 
efficient C in Eq. (Q) is assumed C ~ ^-ioi^Mq)^ ■ The 
mass of a rich cluster is 5-10 14 -5-10 15 Mq, 2-7% of it 



(see, e.g., Voloshchuk 1984). Numerical solution of the 
Smoluchowski equation gives £2 ~ 0.02, £ 4 / 3 w 0.26, 
£3/2 « 0.1 for U = c 2 (M 1 + M 2 ) 2 , U = c 4/3 (Mi + M 2 ) x 
(M* /3 + M 2 1/3 ), U = c 3 / 2 (M 1 + M 2 )(M 1 1/2 + M 2 1/2 ) re- 
spectively (if the initial distribution $0 has a tail, £ u may 
essentially depend on $0 and be much less). Assuming 
No ~ M /Mo where M is an average density of the mass 
contained in galaxies and expressing the variables in as- 
tronomical units, we obtain the order of magnitude for t cr : 

( 2-10 15 v^M~ 1 (M/p)- 1 years (it = 2), 
t CI ~ i 4-10 13 u 7 A/ ( T 1/3 (X/p)- 1 years (u = 4/3), (19) 
I 10 14 VrMfT^iM/p)- 1 years (u = 3/2). 

Here v 7 = tw/(10 7 cm s" 1 ), M 6 = M /(10 6 M ), M/p 
is the ratio of the local density of the mass contained in 



is contained in galaxies (Bohringer 1995). The size (of 
the central part) being of the order of one megaparsec, 
the ratio M/p may be several hundred to several thou- 
sand. Assuming M/p ~ 10 3 we obtain that, for a clus- 
ter with a low velocity dispersion (~ 300 kms _1 ), the 
critical time is less than the age of the Universe on condi- 
tion that masses of the initial galaxies (which then merge) 
M ~ 10 9 -10 10 Mq or more (for the case M > M b , 
u = 3/2). A dose estimate for Mo can be obtained also 
for the regional M -C Mb, u = 2. On the other hand, we 
may consider the formation of massive galaxies and the 
mass function tail only if the initial mass Mo is much less 
than a typical mass of a large galaxy (~ 10 11 Mq). For a 
cluster with a bigger velocity dispersion (1000 km s" 1 or 
more) it is much more difficult to satisfy condition (|l9|): 
Mb is large, so the kernel with u = 2 should be taken; 
t cr is proportional to u 3 ms and can be less than the age of 
the Universe only for very high density (M/p ~ 10 4 ) and 
large enough initial masses (M ~ 10 10 Mq). Thus, pos- 
sible dependence of u rms on time due to cluster evolution, 
its space nonhomogeneity etc. can essentially influence the 
role of mergers, especially on small masses. 

The estimate for M/p given above is based on the as- 
sumption that dark matter belongs to the whole cluster 
rather than to individual galaxies. If dark matter is con- 
centrated in galaxiesQ, the ratio M/p may increase by an 
order which results in the same decrease of t CT (according 
to Eq. ©). 

So, the conditions necessary for the "explosive" process 
of mergers may be realized in many clusters. 

After a cD-galaxy which contain a significant part of 
the total mass has formed in the cluster centre, the dy- 
namics of the cluster is largely determined by attraction 
to this galaxy, and the model considered in this paper 
breaks down. Besides random collisions, spatial inhomo- 
geneity and mass segregation become essential: due to dy- 
namical friction, most massive galaxies gradually gets into 
the centre and are swallowed by the cD-galaxy. However, 
before t CT , when there is no yet cD-galaxy in the cluster, 
galaxy mergers can be considered as random pairwise en- 
counters with probability given by Eq. (^). It is also clear 
that mass segregation at a late stage of cluster evolution 
should be computed together with mass function evolu- 
tion, using a spatially inhomogeneous Smoluchowski type 
kinetic equation, which is a much more complicated prob- 
lem. 



6 The choice of the region depends on the relation between 
Mo and Mb. If Mo Mb then the case u = 2 is realized; if 
M ^ Af b then u = l + /3. 

7 At least cD-gala xies m ay contain a large quantity of dark 
matter (Ikebe et al. 1996| ). 
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Note that galactic mergers due to dynamical friction 
were discussed by Hausman & Ostriker ( 1978) ) . In spite of 
the difference in the model, the algorithm for simulation of 
mergers, used in their work, is analogous to the algorithm 
we use (and, thus, equivalent to the Smoluchowski equa- 
tion), though the expression for the merger probability is 
different (U oc MiM 2 , which also gives the "explosive" 
evolution) . 

As shown in Sects. |^-|| (see also Kontorovich et al. 
1995b} Krivitsky |1995| ), the mass function formed by 
mergers with the probability given by Eq. (m) is rather 
steep (a » 2 for U oc (Mi + M 2 )(Aff + m[), a ~ 2-3 
for U oc (Mi + M2) 2 ; in the latter case the asymptotics 
seem to be non-power, so a cannot be determined accu- 
rately) . It is quite possible that the obtained values of the 
slope correspond to the steepening of the cluster galaxy 
luminosity function at the faint end, which was recently 
discovered (de Propris et al. 1995; Kashikawa et al. 1995} 



Bernstein et al. 1995). According to observational data 



analysed in these works, the effective value of the slope for 
faint galaxies in clusters increases up to 2-2.2 (though, as 
Bernstein et al. note, there is a different interpretation of 
these results) . Possibly, this part of the luminosity func- 
tion is formed by mergers and can be described by the 
intermediate asymptotics for the Smoluchowski equation 
(if the latter can be extended to small enough masses). 
Recent HST data shows also an exce ss of f aint low-mass 
objects for field galaxies (Cowie et al. 1995|) . 

The appearance of relatively steep intermediate 
asymptotics (a » 21 can be easily understood from the 
following arguments!! Both obtained values for the index 
(a w 1.9 for u = 4/3 and a w 2.1 for u — 3/2) are 
within the range (u + 2)/2 to (u + 3)/2. The mass func- 
tion-, with a = (u + 3)/2 corresponds to a constant mass 
fluxu to infinity (i.e., to cD-galaxy, in-our case). However, 
due to nonlocality of the distributionEl (\u2 — Ui\ > 1, see 
Vinokurov & Kats 1980) such a solution is not realized 
exactly in both our cases. Nonlocality leads to an essen- 
tial role of interactions between low-mass and high-mass 
galaxies. Then the number of massive galaxies is approx- 
imately conserved, and the constant flux of their number 
to in finity corresponds to (u + 2)/2 (Kontorovich et al. 
1993 ). Since none of these limit cases is realized, the index 
is situated between these values: 1.67 < a w 1.9 < 2.17 
(u = 4/3), 1.75 < a « 2.1 < 2.25 (u = 3/2), and is 
rather close to their arithmetic mean value (as we can see 
both from the simulation and the numerical solution of 
the Smoluchowski equation). 



We consider here the case U oc (Mi+M 2 )(Mf +M%), when 
the mass function is close to a power law. 

9 Solutions with a constant flux of a conserved variable are 
analogous to Kolmogor qv spe ctra in the weak turbulence the- 
ory, see Zakharov et al. 1992| . 

Nonlocality corresponds to divergence of the collision inte- 
gral for the power-law distribution. 



The density ratio M/ p in the above estimates was one 
of the most important parameters which control the pos- 
sibility of effective merger process. The local value of this 
parameter may vary in a very wide range: from 1 (scales 
exceeding an average distance between massive field galax- 
ies) to 10 7 (if we take an average density of a galaxy 
~ 10 -22 g cm~ 3 for M.). As was shown above, in clusters 
this parameter is large enough (~ 10 3 or even 10 4 ) to yield 
an "explosive" evolution due to mergers. Local concentra- 
tions may enable analogous phenomen a for field galaxies 
at large z (see, e.g., Komberg & Lukash 1994; Kontorovich 



1994) 



Morphological changes in cluster galaxies, which is one 
of the results of mergers, may be related to the chang e 
of the angular momentum distribution (cf. Toomre 1977[ ) . 
The possibility of dependence between Hubble's morpho- 
logical type and an effective angular momentum has been 
discussed in the literature (see Fig. 1 in Polyachenko et al. 



1971) and confirmed by an independent analysis of ob- 
servational data (Kontorovich & Khodjachikh 1993; Kon- 
torovich et al. 1995a). However, this dependence needs 



special consideration which is beyond the scope of this 
work. 

The above consideration of cluster evolution takes into 
account only galaxy mergers in a spatially homogeneous 
model. It allows to obtain the "explosive" evolution, the 
steep part of the mass function, cD-galaxies, rapid evolu- 
tion of galaxy morphological types, and a mean value of 
the dimensionless angular momentum which does not de- 
pend on the details of the initial distribution. In the same 
time, this approach has obvious limitations. The "explo- 
sive" evolution does not produce Schechter's mass func- 
tion with a w 1.25. It is possible that the effective merger 
probability U changes at a late stage of cluster evolution 
(when massive galaxies have formed) in such a way that 
u c ff becomes, less than one and the "explosive" process 
slows downliil, which leads to a flatter $(M). Another pos- 
sibility is that this part of the mass function may not be 
formed only by mergers. 
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